QCD Equations of State and the QGP Liquid Model 
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Recent advances in the study of equations of state of thermal lattice Quantum Chromodynam- 
ics obtained at non-zero baryon density allow validation of the quark-gluon plasma (QGP) liquid 
model equations of state (EoS). We study here the properties of the QGP-EoS near to the phase 
transformation boundary at finite baryon density and show a close agreement with the lattice results. 



PACS numbers: 12.38.Mh 

The ab-initio exploration of Quantum Chromo dynam- 
ics (QCD) on the lattice at finite temperature and baryon 
density has seen rapid recent progress 0, 0, 0- It is 
already known that for vanishing baryon density lat- 
tice results can be described in terms of ideal quantum 
gases allowing only for lowest order perturbative inter- 
actions, provided that a non-perturbative temperature 
dependence of the couplingconstant is introduced, along 
with vacuum latent heat We will refine this quark- 
gluon plasma model, and demonstrate that it agrees ex- 
tremely well with the finite baryon density lattice results, 
except in a region of temperature in vicinity of the phase 
transition/transformation domain. 

A strong motivation to have an understanding of the 
thermal QCD lattice results in terms of a set of more in- 
tuitive degrees of freedom originates in the need to treat 
a fast evolving system created in relativistic heavy ion 
collisions. At the onset of the QCD thermal matter for- 
mation, e.g., at To — 0.5 fm, one cannot expect a full 
chemical equilibrium to have been reached, i.e., the phase 
space occupancy of quark and/or gluon gases has not ap- 
proached unity. In a model as we will present these situa- 
tions can be easily incorporated, allowing study of initial 
temperature, phase space occupancy evolution, and other 
physical conditions in the realistic environment of high 
energy relativistic heavy ion collisions. 

More generally, a handy model of QCD matter equa- 
tions of state employing established physical concepts al- 
lows to study many physical properties which are only 
difficult to infer from the numerical lattice results. But 
perhaps the most important question we pursue is if 
one can already in the temperature domain explored by 
present day experiments describe thermal QCD matter 
in terms of nearly free quark and gluon degrees of free- 
dom. If this is confirmed, we would be able to verify 
if the physical system formed in these reactions is the 
deconfined quark-gluon plasma. 

A systematic perturbative expansion within the frame- 
work of thermal field theory of the interacting quark- 
gluon gas converges poorly p|. The situation is illus- 
trated for the case of a pure gauge SU(3) case in figure^ 
where we see that the ratio of the computed perturba- 



tive pressure P to the Stefan-Boltzmann pressure Pq is 
oscillating around unity depending on the order n of the 
expansion in the coupling constant g considered. Con- 
vergence seems to occur only in the asymptotic freedom 
limit of relatively high temperature T, beyond the exper- 
imental reach. Better agreement with the lattice results 
(solid dots) can be achieved (solid line) at a relatively low 
T, here expressing the temperature in units of the renor- 
malization scale A^jg, provided that the unknown relative 
coefficient of the g term is optimized to the value 0.7 . 

We advocate a very different approach, combining the 
perturbative result with key nonperturbative features. 
Our QGP liquid model arises from the empirical observa- 
tion that the lowest order thermal perturbation contribu- 
tion evaluation of the QCD matter properties, combined 
with a non-perturbative temperature dependent strong 
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FIG. 1: Ratio of the pure gauge pressure with Stefan- 
Boltzmann pressure P/Po as function of temperature T in 
units of renormalization scale A^jg-, results of Ref. [fj. 
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coupling constant agrees with the QCD thermal lattice 
results, once the other often used non-perturbative effect, 
e.g. bag constant, has been incorporated. 

The QGP-liquid partition function is assumed to have 
the form, 
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and: 



c 3 = 1 - — . (2) 
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e^ /T . A value B 



We recall that fi b = 3fi q and X q 

(0.211 GeV) 4 fits best the lattice results we consider here. 

The temperature dependence ct s (T) is obtained from 
a s (x) setting the energy scale x to be: 



27r/r 1 WlT^T? A q = 2^(7r/ 



(3) 



a s (x) is obtained integrating the renormalization group 
equation, incorporating physical thresholds for heavy fla- 
vor. Use of semi-analytical formulas with fixed active 
number of quarks introduces an unacceptable error in 
the value of a s and has lead to false conclusions in some 
earlier work. We evaluate: 
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02 ert is the beta-function of the renormalization group in 
two loop approximation, and 
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2^ 



bi = 



51 - 19ra f /3 
4^ 



/32 Mt does not depend on the renormalization scheme, 
and solutions of Eq. Q differ from higher order renormal- 
ization scheme dependent results by less than the error 
introduced by the experimental uncertainty in the mea- 
sured value of a s (/i = M z ) = 0.118 + 0.001-0.0016, used 
as the initial value in numerical integration of Eq. Q . We 
note for // q — > the empirical form: 
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T c = 172MeV. (5) 



1 + 0.721n(T/T c ) 

When we explore the effect of the finite quark masses, 
we introduce the correction S m In Zq GP to the partition 
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where i = q,s. We use m s = 170 MeV. Note that 
we did not allow for the QCD-thermal a s correction in 
5 m In Zqqp- When we compare withe lattice results, we 
use the values m q = 65 MeV and m s = 135 MeV, as re- 
ported to have been used in lattice simulations we com- 
pare with. Strange quarks enter below only when we 
consider the absolute pressure. Similarly, when we ex- 
plore how an effective gluon mass alters the pressure, we 
consider a correction 
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We need to introduce a finite thermal glue mass rriQ ~ 
0.2 GeV in order to reach for T > T c a line width agree- 
ment with lattice results. This is not the thermal gluon 
mass which, as we have discussed elsewhere, can be used 
as an alternative way to express the perturbative QCD 
effect p. 

In order to be able to consider our results in currently 
ongoing experiments, it is important that the QGP-liquid 
equations of state are verified at chemical potential which 
is relatively small, up to /!;, ~ 100 MeV at initial condi- 
tions. However, wc extend this study to all values consid- 
ered in the lattice simulation [lj, fib = 100,210,330,410 
and 530 MeV, which aside of RHIC, also encompasses 
the physical reach of CERN-SPS experiments. Our ob- 
jective is to test if for T > 1.5T C the model and lattice- 
'experiment' agree, while near to the critical tempera- 
ture domain significant modifications due to the non- 
perturbative features of QCD must be expected. 

Physical properties we consider will be rendered di- 
mensionless by considering a suitable ratio with cither 
T 4 (pressure P) or T 3 (baryon density To asses if 
the liquid model has the right density dependence we 
consider in turn the pressure P modification at finite 
baryon density, AP = P(T,fi b ) - P(T,fi b = 0), the 
baryon density n b , the pressure P(T,fi b = 0) and finally 
e(T, fib) — 3P(T, fit), which vanishes for ideal gases. In the 
following figures all 'experimental' points are taken di- 
rectly from the postscript format of the figures of Ref. , 
into which we add our results. 

The reader should note that for T < T c equations of 
state with hadrons must be employed. Hence we do not 
here include results for this domain of lattice results. The 
study of this confined phase domain involves also mod- 
eling of excluded volume effect, and/or within the boot- 
strap model, the understanding of the singularity of the 
hadron mass spectrum. Given the availability of lattice 
results we hope to return to discuss these intricate mat- 
ters in near future. 

The change on AP at finite baryochemical potential 
is shown in figure [3 We note that the agreement with 
the lattice results at large T is very satisfactory. This 
result depends on the empirical choice made for the de- 
pendence on chemical potential of the scale of a s , see 
Eq. Had we omitted fi q , or doubled the coefficient 
of fi q , there would be a well visible deviation in figure 



3 



0.6 



0.2 



[ 1 1 111 1 1 1 [ 

\ 


I I I | 1 1 


1 1 | I I 


i i 










i \ 
















■ *\ \ 






- 
































i 'Sn V 








■ y 

■ / « * 









1 1.5 2 2.5 3 

T/T c 

FIG. 2: AP ee P{T,fi b ) - P(T,/x 6 = 0) normalized by T 4 
as function of T/T c for ft 6 = 100, 210, 330, 410 and 530 
MeV from bottom to top. Data points from ref. .1], solid 
lines massless liquid of quarks. Dashed (and mostly invisible) 
results with finite mass correction applied for m q — 65 MeV 
as used in lattice data. 



12 Thus the empirical choice of the Matsubara frequency 
as the combination of temperature and chemical poten- 
tial is qualitatively verified by lattice results. For small 
chemical potential ^ = 100 MeV, we see that there is 
agreement for T > 1.1T C , at a baryochemical potential 
Hb = 530 MeV, agreement with lattice data is assured for 
T > 1.5T C . The practically invisible dashed lines (hidden 
mostly under solid lines) show the here negligible effect 
of finite quark masses. By definition, AP depends only 
on light quark degrees of freedom. The agreement we see 
is thus not testing strange quark, or glue behavior. 

A similarly remarkable agreement is obtained for the 
baryon density, shown in figure [3J Here a finite quark 
mass used in lattice simulations becomes visible (dashed 
lines). Once we allow for it, we see agreement within 
the lattice data error for T > 1.2T C for all baryochemical 
potentials. For very small baryon density we can expect 
agreement down to near the critical temperature. 

Since both the pressure difference AP and baryon den- 
sity Ub show satisfactory behavior as function of chemical 
potential, it is expected that the change of energy density, 
and entropy, with chemical potential is well described, 
both being derived of AP with respect to j3 — 1/T and 
T and suitable linear combination with ub- However, 
in order to fully certify the liquid model, we need now 
to fine-tune the behavior of pressure P at zero chemical 
potential. 

In figure 0J we show the pressure P(T, fij, — 0)/T 4 
as function of temperature T/T c . We find that there is 
some difference in the lattice result shown in and ear- 
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FIG. 3: Baryon density ub normalized by T 3 as function of 
T/T c for fi b = 100, 210, 330, 410 and 530 MeV from bottom 
to top. Data points from ref. Q], solid lines massless liquid of 
quarks. Dashed lines: allowance is made for m q = 65 MeV as 
is used to obtain the lattice data. 

lier work , which practically agreed with our solid line 
result in figure 01 To obtain the QGP model agreement 
with the more recent lattice data [J we need to: 
a) change the value of the Bag constant, from our earlier 
value (0.195 GeV) 4 to B = (0.211 GeV) 4 , in response to 




FIG. 4: The pressure P(T, ^ = 0) normalized by T 4 as func- 
tion of T/T c . Data points from ref. 0. Solid lines: massless 
gluons with B = (0.211 GeV) 4 . Dashed line allows for a finite 
mass m G = 200 MeV. 
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FIG. 5: (e - 3P) /T 4 as function of T/T c for /i b = 0, 210, 410 
and 530 MeV, for B = (0.211 GeV) 4 and T c = 173 MeV. Data 
points from ref. Solid lines: massless gluons. Dashed 

lines: allowance is made for ma = 200 MeV. All chemical 
potential lines coincide within line-width. 

somewhat higher value of T c , see solid line in figure 01 
b) introduce an effective gluon mass ma — 0.2 GeV in 
order to reduce the number of effective degrees of free- 
dom as is shown by the dashed line in figure 0] 
This procedure does indeed produce the expected agree- 
ment within the line-width. 

We do not have a lattice entropy figure to compare 
with, but we are assured of a valid result by our abil- 
ity to reproduce the shape of the pressure functions, see 
figures |3 and 0] A sensitive test of this assertion is ob- 
tained considering (e — 3P)/T 4 as function of T/T c in 
figure [SJ All finite chemical curves we plot coincide and 
hence only one is visible in the diagram. This agreement 
with lattice results for T > 1.1 5T C confirms that we have 



obtained a remarkably precise representation of the be- 
havior of equations of state of deconfined QCD matter, 
except in direct vicinity of the critical temperature. 

We have shown here that QCD equations of state at 
nonzero baryon density for T > 1.5T C , but even very near 
to T = T c for small chemical potentials, behave in a way 
expected for the quark-gluon plasma phase. The remark- 
able agreement has been made possible by introduction of 
the first order thermal interactions which comprise pre- 
cise nonperturbative thermal coupling a s (T, (i q ) and the 
vacuum latent heat B. In the glue sector, effective gluon 
mass me — 200 MeV is also required to obtain exactly 
the expected number of degrees of freedom, required to 
model the absolute magnitude of the pressure to within 
line-width. It is important to realize that even without 
such detailed refinements which may indeed change again 
as the lattice results evolve, the approach we advocate 
does not suffer from the convergence problem illustrated 
in figure ^ an d it provides a very good and natural de- 
scription of the wealth of the lattice results. 

We conclude that thermal lattice QCD matter behaves 
just like quark gluon plasma, and thus a naive use of 
the QGP model is appropriate, provided that suitable 
coupling strength a a {T,fi q ) and vacuum bag constant 
B is introduced. Moreover, considering that at RHIC 
the following initial conditions have been reached 
T > 1.5T C ~ 260 MeV for \ q = 1.09, corresponding to 
the initial baryochemical potential [ib = 3T In \ q ~ 45 
MeV, the results presented confirm that the QGP state 
is established in these reactions. Moreover, we are re- 
assured that we can explore in detail the RHIC initial 
conditions, allowing for chemical composition dynamics. 
This should lead us to understanding of QGP properties 
and conditions in RHIC reactions. 
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